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Sputtering in oblique C-type shocks 
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ABSTRACT 

We present the first results for the sputtering of grain mantles and cores obtained with self- 
consistent multifluid hydromagnetic models of C-type shocks propagating through dusty me- 
dia. The threshold shock speed for mantle sputtering is about 10 km s _1 and is independent 
of density. The mantles are completely vapourised in shocks with speeds of 20-25 km s _1 . 
At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruc- 
tion is not total in any C-type shock because grains are not completely destroyed in shocks 
with speeds near the minimum speeds at which J-type shocks appear. Due to the density- 
dependence of the critical shock speed for this transition, higher gas-phase SiO fractional 
abundances are produced behind shocks propagating in lower density gas. For shock speeds 
near the threshold speeds for both core and mantle sputtering, sputtering is much greater for 
shock velocities at smaller angles relative to the upstream magnetic field. At higher shock 
speeds, the angular variation is still present but less pronounced. 
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1 INTRODUCTION 

A number of computational studies of sputtering in self- 
consistent models of perpendicular C-type hydromagnetic 
shocks in dusty, dark regions have been performed (e. g. 



Flower & Pineau des Forets 


19951; Schilke et al.ll 19971: iMav et al.l 


200d:lGusdorf et alj|2008bl: 


Guillet et al. 201 1 ). However, the onlv 



numerical investigation of sputtering in oblique C-type shocks 
in such regio ns is not based on self-consistent shock models 



dCaselli et al.lll997h. 



Pilipp & HartquisJ J 1994) attempted to construct a self- 



consistent, four fluid model of an oblique C-type shock in which 
dust grains are described as a charg ed fluid, as they are in the per- 
pendicular C-type shock mod els of Pilipp et all 1 1990l). However, 
the integrations performed by IPilipp & Hartquist J 1994) failed to 
yield solutions corresponding to steady fast-mode shocks. Iwardld 
i 19981) pointed out that the failure is due to the saddle -point nature 
of the downstream state of a fast-mode dusty shock. iFallej d2003T) 
developed a numerical scheme for solving the time-dependent mul- 
tifluid hydromagnetic equations that enables the construction of 
steady fast-mode solutions even when local equilibrium does not 
obtain. IVan Loo et al. have used the scheme along with non- 

equilibrium chemistry and treatments of grain charge and heating 
and radiative losses to investigate the evolution of oblique C-type 
shocks propagating in to media with constant upstream properties. 
lAshmore et alj feOld) extended the applications to oblique C-type 
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shocks propagating through media with varying upstream proper- 
ties. 

Here we use models like those of IVan Loo et al .1 ( 120091) to 
study sputtering in oblique C-type shocks in dusty, dark molecu- 
lar regions. Section 2 contains descriptions of the assumptions, the 
physical processes included, the chemical network, the evaluation 
of the sputtering rates and upstream boundary conditions. Results 
are presented in Section 3, and Section 4 contains a discussion of 
the results and conclusions. 



2 THE MODEL 

The model includes neutral, ion, electron and two grain fluids. 
Flow variables depend on only the spatial coordinate x and time 
t. The upstream magnetic field has components in the x and y di- 
rectio ns, and the shock prop agates in the x direction. Equati ons 1 
- 8 of Ivan Loo et al] J2009T ) are solved with the scheme of iFallel 
J2003h . The inertia of the charged fluids, including the grain flu- 
ids, is neglected in this approach, so that their dynamics is solely 
determined by the balance between the Lorentz force and collision 
forces with neutral particles. Additionally, the grain fluids are as- 
sumed to have zero pressure. The average charge on a grain is de- 
termined by the condition that the magnitude of the current of ions 
onto a grain equals the magnitude of the curr ent of electrons onto a 
gram; Eqs. 10 and 11 of lHavnes et "all i 19871, see also Appendix lAl> 
are us ed to calculate these currents. The numerical method of lFalld 
d2003h provides a rigoro us means to model stea dy shock struc- 
tures and, as discussed by (Van Loo et al ,l2009h and lAshmore et all 
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Figure 1. The fractional abundance of H2O produced by gas-phase chem- 
istry (solid) and mantle sputtering (dashed) as function of the shock speed 
for rajjo = 10 6 cm~ 3 . 



( |201(]|) , reproduces the results of, among others, Draine et"al] 
]l983n.|Pilipp & Hartquistl h994h.lcMpman & Wardld l l2006l) and 
iGuillet. Pineau des Forets. & Jones! J2007h . Furthermore, although 
it neglects grain ine rtia, the me t hod can accurately m odel neutral 
subshock structures dFallel2003l : lAshmore et al.ll2010h . 

The evolution of the number density, n(X), of chemical 
species X is governed by a rate equation of the form: 



dn(X) + dn{X)u s = 



Of 



dx 



(1) 



u 3 is the x component of the ion or neutral velocity depending on 
whether X is an ion or a neutral. s(X) is the net rate of forma- 
tion or destruction of species X per unit volume per unit time due 
to the gas phase reactions and grain surface recombinations listed 
in Table Q] and the return of material to the gas phase due to sput- 
tering. Though no reactions involving SiO are shown in Table Q] 
atomic Si returned to the gas phase due to sputtering is assumed to 
be converted to SiO immediately. 

The model includes the radiative losses due to O, CO, H2 
and H2O, which are calculated for the abundances of O and 
H2O given by the integration of the chemical rate equations 
and the assumed constant fractional abundances of CO and H2. 
Given th ose abundance s, the c ooling rates are evaluated as they 
were by Van Loo et al.l d2009l) . They adopted the treatments of 
iDraine et alj 1 19831) for the cooling of ele ctrons and for 01 emis- 
sion, used data from lHartquist et all dl98Clh for the cooling of neu- 
trals by H2 and implemented the treatment of CO cooling due to 
iHollenbach & McKeel d 19791) and that of H2O cooling developed 
bv lNeufeld & Melnicldll987l) . 

Sputtering due to H2, CO, O, Mg, H 2 and SiO projectiles 
is included. As mentioned above, the return of material to the gas 
phase due to sputtering contributes to the source terms in the rate 
equations governing the gas phase chemistry; these contributions to 
the source terms are particularly important for the model gas phase 
abundances of Mg, H2O and SiO. A grain consists of a refractory 
core made of olivine (MgFeSi04) and a mantle of water ice. A core 
remains intact until all of the ice has b een sputtered. The sputterin g 
rates are calculated as they were by Jimenez-Se rra et al. I J2008h . 
who used Eq. 27 of IDraine & Salpeta 1 19791) for the sputtering 
rate. F or all projectiles, they used Eq. 31 of IDraine & Salpeterl 
dl979h to obtain the sputtering yields of ice at normal incidence. 
They used the same equation to obtain the corresponding yield 



for core material colliding with H2 . Ho wever, for projecti les other 
than H2, they employed the results of iMav et al. I feood) for the 
sputtering yields of core material at normal incidenc e. Following 
IDraine & Salpeterl dl979l) . |jimenez-Serra eTal] d2008l) took the an- 
gle averaged value of a yield to be twice the value of the normal 
incidence yield. 

The upstream magnetic field strength is Bo = 
l/iG(nHo/l cm -3 ) 1 / 2 , where nHo is the upstream value of 
nH, the number density of hydrogen nuclei. A shock propagates 
with a speed v s at an angle 6 with respect to the upstream magnetic 
field (for an oblique shock, 0° ^ 6 <c 90°). The grains have 
radii of 0.4 /j,m and 0.04 fim and masses of 8.03 xl0~ 13 g and 
8.03 xl0~ g. In the upstream region, one percent of the mass 
is contained in grains with the number density of smaller grains 
316 times that of the l arger grains ( i .e. th e grains follow the 
grain-size distribution of iMathis et alj dl977l) ). As all the grains 
are silicate-based in our model, the fractional abundance of Si, 
i.e. 5.8 x 10 ~ 5 , is somewhat higher than the cosmic value of 
3.37 x 10~ 5 dAnders & Grevessdll98^) . The fractional abundance 
of CO to hydrogen nuclei has a constant value of 5 x 10~ 5 
ahead, throughout and behind the shock, and the upstream gas 
phase fractional abundances of O, Mg, H2O and SiO are 4.25 
x 10~ 4 , 10 -7 , and 0, respectively. Note that, including the 
fractional abundance of CO and of olivine, the total abundance 
of O in our model exceeds the cosmic abundance of O by less 
than a factor of two (the lates t measure d value being 3.2 x 10 -4 
dMever. Jura. & Cardelli|[l998l) although Isiiow & Witj Jl996T) give 
4.6 x 1 0~ 4 ). The init i al frac tional ionisation is calculated as it 
was by I Van Loo et all d2009l) : HCO + and Mg+ are initially the 
only ions, but the chemistry leads to the production of significant 
fractional abundances of and H,-iO + before heating by the 
shock occurs. Although upstream of a shock there is no water 
in the gas phase, the upstream fractional abund ance of water ice 
molecules on the grains is 7.25 x 10" 5 (cf. IWhittet & DuTevl 
1 199 ll) : the upstream value of the ratio of the total mass of ice on 
small grains to the total mass of ice on big grains is taken to be 
3.16. 

We followed lVan Loo et al. by assuming that the initial 

flow corresponds to that of a J-type shock. The presented results are 
for the steady flows yielded by time integrations. 



3 RESULTS 

We focus on the abundances of H2O and SiO as mantle and core 
sputtering releases these molecules into the gas phase. However, 
water is also produced by shock chemistry. To quantify the contri- 
butions of both processes we discuss the chemistry first and then 
grain sputtering. 



3.1 Shock chemistry 

The specified initial relative abundances are not in steady state 
chemical equilibrium. The upstream gas phase ion abundances ini- 
tially evolve on a timescale short compared to the flow time through 
the shock and then evolve further as the abundances of other key 
species (e. g. H2O) vary. HsO + and HCO + are the most com- 
mon molecular ions, but the Mg + fractional abundance is about 
an order of magnitude greater than the total fractional abundance 
of molecular ions. The ratio of the total atomic ion abundance to 
the total molecular abundance depends on the assumed elemental 
abundance of metals. 
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Table 1, The chemical network. The column titled "Num." lists the reaction number, that labelled "Reac- 
tion" gives the reaction partners and products for each, that designated "Rate or Rate Coefficient" lists the 
a dopted rates and ra te coefficients and that labelled "Ref." gives the source of each coeffi cient: Source 1 
is lPilipp et ai] ll990l) . 2 is the UDFA database described bv lWoodall eTal] fe007l) . 3 is lGuillet et all 120091) 
and 4 indicates that the value is calculated with a routine in the code. 



Num. 


Reaction 


Rate (s 1 ) or Rate Coefficient (cm 3 s 1 ) 


Ref. 


1 a 


r±2 + CR — > + e + CR 


17 

3 x 10 


1 


2 


H]}~ + CO -> HCO+ + H 2 


1.7 x 10 -9 


1, 2, 3 


3 


Hg" + H2O — > H3O+ + H2 


5.9 x 10 -9 


2, 3 


4 


Hg" + grain - — y H2 + H + grain 


(Ti„ + r,-,) n(Hi)/nj 


4 


5 


+ -> H3O+ 


8.4 x 10 -10 


2. 3 


6 


HCO+ + e -> H + CO 


2.76 x 10 -7 (T e /300) -0 - 64 


3 


7 


HCO+ + Mg -> H + CO + Mg+ 


2.9 x 10 -9 


2 


8 


HCO+ + H 2 ->• H3O+ + CO 


2.5 x 10 -9 


2, 3 


9 


HCO+ + grain" — > H + CO + grain 


(r, 9 + n(HCO+)/n, 


4 


10 


Mg+ + e — > Mg + hu 


2.8 x 10 -12 (T e /300) -0 - 86 


1,2 


1 1 


Mg+ + grain - — > Mg + grain 


(r i9 + T is ) n(Mg+)/n; 


4 


12 6 


Mg+ + H 2 -> MgH+ + H 


3.0 X 10 -10 exp(-27854AT eff ) 


1 


13 


MgH+ + e -> Mg + H 


1.1 x 10 -7 (T e /300) -0 5 


1 


14 


MgH+ + grain - — > Mg + H + grain 


(r i9 + r is ) n(M g H+)/ni 


4 


L5 C 


H3O+ + e H 2 + H 


3.29 x 10 -7 (T e /300) -0 ' 5 


2. 3 


16 


H3O+ + Mg -> H2O + Mg+ + H 


1.8 x 10 -9 


1 


17 


H3O 4 " + grain - — > H 2 + H + grain 


(r is +r is )n(H 3 0+)/ni 


4 


I8 d 


O + H 2 -> H2O 


3.14 x 10 -13 (T n /300) 2 ' 7 exp(-3150/T n ) 


2, 3 



a Reaction 1 is an amalgamation of two reactions: H2 + CR — > H J + e + CR and H^" + H2 — > H^" + H 
where CR denotes a cosmic ray. The adopted coefficient is that of the first reaction since the proton 
transfer reaction occurs sufficiently rapidly. 

b T ef f in the rate coefficient of Reaction 12 defines a weighted kinetic temperature at the ion-neutral drift 

speed and is given by Eq. 43 of lFlower et alj dl985D . 

c Reaction 15 is an amalgamation of: 

H3O+ + e^OH + H 2 (6.02 x 10 -8 (T e /300) -0 - 5 ), 

H3O+ +e ->OH + H + H (2.58 X 10 -7 (T e /300) -0 - 5 ) and 

H3O+ + e ^H 2 + H (1.08 x 10 -8 (T e /300) -0 5 ). 

OH is not treated as a separate species, because the OH produced by the first two reactions quickly reacts 
with H2 to yield H2O (see d ). 

d Reaction 18 is also an amalgamation of two reactions: O + H2 — > OH + H and OH + H2 — > H2O + H. 
The coefficient adopted is that of the first (slower) of the two. 



Table 2. Fraction of mantle that is sputtered and gas phase water fractional abundance enhancement due to sputtering. Results are given for 
large and small grains and for both combined. Although some small grain sputtering does occur in the 9 km s -1 , njjo = 10 6 cm -3 case, the 
fraction of ice removed was < 10 — 4 . 







Lar 


ge grains 


Small grains 


Large 


& Small grains 


Shock Speed 


nHO 


Fraction of 


Contribution to 


Fraction of 


Contribution to 


Fraction of 


Contribution to 


(kms -1 ) 


(cm -3 ) 


ice removed 


downstream H2O 


ice removed 


downstream H2O 


ice removed 


downstream H2O 


10 


1.00(4) 


1.20(-3) 


3.09(-8) 


3.80(-3) 


2.09(-7) 


3.3K-3) 


2.40(-7) 


12 


1.00(4) 


4.00(-2) 


1.04(-6) 


5.50(-2) 


3.03(-6) 


5.61(-2) 


4.07(-6) 


14 


1.00(4) 


1.75(-1) 


3.05(-6) 


2.00(-l) 


1.10(-5) 


1.94(-1) 


1.41(-5) 


16 


1.00(4) 


4.50(-l) 


7.83(-6) 


4.80(-l) 


2.64(-5) 


4.72(-l) 


3.42(-5) 


20 


1.00(4) 


1.00 


1.74(-5) 


1.00 


5.51(-5) 


1.00 


7.25(-5) 


8 


1.00(6) 


0.00 




0.00 




0.0 




9 


1.00(6) 


0.00 




<1.00(-4) 


« 0.00 


< 1.00(-4) 


Ri 0.00 


10 


1.00(6) 


0.00 




5.00(-3) 


2.75(-7) 


3.79(-3) 


2.75(-7) 


12 


1.00(6) 


2.50(-3) 


4.36(-8) 


6.95(-2) 


3.83(-6) 


5.34(-2) 


3.87(-6) 


14 


1.00(6) 


2.50(-2) 


4.36(-7) 


2.63(-l) 


1.45(-5) 


2.06(-l) 


1.49(-5) 


15 


1.00(6) 


5.00(-2) 


8.71 (-7) 


4.25(-l) 


2.34(-5) 


3.35(-l) 


2.43(-5) 


16 


1.00(6) 


7.78(-2) 


1.36(-6) 


6.07(-l) 


3.34(-5) 


4.79(-l) 


3.48(-5) 


18 


1.00(6) 


1.8(-1) 


3.14(-6) 


1.00 


5.51(-5) 


8.03(-l) 


5.82(-5) 


20 


1.00(6) 


3.4(-l) 


5.93(-6) 


1.00 


5.51(-5) 


8.42(-l) 


6.10(-5) 


25 


1.00(6) 


1.00 


1.74(-5) 


1.00 


5.51(-5) 


1.00 


7.25(-5) 
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Figure 2. Fraction of ice ejected from grain mantles for a shock speed of 16 km s _1 , nno = 10 4 cm -3 and 8 of 30° (left), 45° (middle) and 60° (right). The 
top row gives the ion (black, dotted), electron (red, dashed) and neutral (blue, solid) temperatures, as well as the neutral density normalised to its upstream 
value (green, dash-dotted). The middle row shows the fraction of ice ejected from large (blue, dashed) and small (red, dotted) grains. The bottom row displays 
the speed of small (red, dotted) and large (blue, dashed) grains relative to neutrals in km s — 1 . 



The assumed upstream water fractional abundance is about 
three orders of magnitude below the water fractional abundances 
obtai ned with standard i o n-molecule chemical models, i.e. ~ 10~ 6 
(e.g. iMillar et alj 1 1991 iNeufeld et all 1 19951) . However, the up- 
stream value is close to the water abundances of ~ 10 -8 inferred 
for quiescent mo l ecular gas towards Orion, Ml 7 and other sources 
by iBergin et al l J200fjv Recent Herschel observations of ortho- 
H 2 in a dark cloud dCaselli et alfcOldh imply values of ~ 1(T 10 . 

Within the shock front the gas heats up, and some reactions 
with barriers and endothermic reactions occur. Also, the heating 
of electrons leads to enhanced electron impact rates on the grains, 
and the drift of the ions relative to the grains contributes further to 
efficient ion recombination on grains, reducing the fractional ion- 
isation within the shock. As usual, the primary sequence for H2O 
formation in the shock involves the reactions of O and OH with 
H2. Figure[T]shows the the fractional abundance of water produced 
in the shock as a function of the shock speed. At the lowest shock 
speed (i.e. 5 km s™ 1 ), the shocked gas does not reach high enough 
temperatures to cause the conversion of much oxygen to water. At a 
shock speed of 10 km s _1 a few percent of the gas phase oxygen is 
converted to water, but it is only at shock speeds close to 20 km s - 
and above at which the atomic oxygen abundance falls by several 



orders of magnitude, and nearly all of the gas phase oxygen is in 
water. The H2O abundance is then a few times 10~ 4 . 

3.2 Mantle sputtering 

Within a C-type shock, grains experience high-speed impacts with 
gas particles. Sputtering leads to a continuous erosion of the icy 
grain mantle by which water is returned to the gas phase. 

Grain mantle sputtering occurs for impact speeds > 
1Q/-JJL km s _1 with /1 t he mean mass of the impacting species 
in atomic mass units (e.g. iMav et alJ|2QOCh . Since CO is the most 
massive abundant projectile with fj, = 28, mantle sputtering starts 
for grain-neutral drift speeds as low as ~ 2 km s . (Note from 
Fig.[2]that appreciable sputtering only occurs at drift speeds above 
6 km s _1 ). This value of the drift speed is independent of the grain 
size, because the mantle molecules are bound with the same energy. 
However, the threshold drift speed is reached at different shock 
speeds for the small and large grains. Small grains have Hall pa- 
rameter^] that are large enough that they move with the ions and 
electrons. The Hall parameter for the large grains is above unity but 

1 The Hall parameter is the ratio between the gyrofrequency of a charged 
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Figure 3. The fraction of mantle ice (blue) and core silicon (red) released 
in the gas-phase by sputtering as function of shock speed for values of njjo 
of 10 4 cm~ 3 (solid) and 10 6 cm" 3 (dashed). The squares show the sput- 
tering fractions for shock angles of 30, 45 and 60 degrees for selected shock 
speeds; the sputteringefficiency decreases with increasing angle. 



sufficiently small that they are only moderately coupled to the mag- 
netic field and, thus, move at velocities between those of the elec- 
trons and the neutrals. Consequently, the grain-neutral drift speed 
for the large grains is smaller than for the small grains at a given 
shock speed. 

No projectile induces noticeable sputtering when v 3 = 5 km 
s - , but when v s has a higher value, small grains are sputtered. 
Table [2] as well as Fig. [5] displays the fraction of the initial man- 
tle ice that is sputtered by a shock for which 9 = 45°. Results 
are given for various values of the shock speed and values of nHo 
of 10 4 cm~ 3 and 10 6 cm -3 . It also gives the enhancement in the 
downstream gas phase water fractional abundance due to sputter- 
ing. If we define the sputtering threshold shoc k speed as the shoc k 
speed for which 1 % of the ice is removed (cf. lCaselli etahll 19971) , 
the threshold shock speed for small grain sputtering is 10 km s _1 
and for large grain sputtering it is 10 (12) km s for an upstream 
number density of 10 4 (10 6 ) cm . 

One also sees that the fraction of material sputtered from small 
grains is not very sensitive to nHo, but particularly at lower shock 
speeds the fraction of material sputtered from the large grains does 
depend on nHo • This again is a result of the coupling of the grains 
to the magnetic field (see above). The large grains have a larger 
Hall parameter in the lower density regime and, thus, have a larger 
neutral-grain drift speed than in the high-density regime. 

For shocks for which 9 — 45°, the highest sputtering rate per 
unit time per unit volume of mantles is due to collisions with CO for 
shock speeds less than about 10 km s _1 and to collisions with H2O 
for speeds above about 14 km s _1 . At intermediate speeds, colli- 
sions with O also contribute non-negligibly as the relative impor- 
tance of CO and H2O collisions varies. Above 30 km s _1 , H2 dom- 
inates the sputtering of both grain types. However, at these shock 
speeds, the icy mantles of both grain species are completely eroded. 
The contribution of the grain sputtering to the water abundance is 
even then negligible compared to the contribution of the gas-phase 
oxygen conversion (see Fig.[TJl. 



particle and the collisional momentum transfer frequency of the charged 
particle with neutrals. 



3.3 Core sputtering 

Once the grain mantles are vapourised, the grain cores begin to 
erode as long as the impact speeds of the projectiles are high 
enough. As the binding energy of a particle in the grain core is 
higher than in the mantle, higher grain-neutral drift speeds are nec- 
essary. Again CO is the sputtering projectile of most importance 
at lower speeds, and small grains start to contribute to the SiO 
gas-phase abundance at shock speeds of 25 km s -1 (see Fig. [3] 
fo r shocks for which 9 — 45°). This result is similar to the one 
of ICaselli et ail i 19971) and iMav et ail d2000h . Collisions with CO 
dominate the sputtering at shock speeds up to 30 km s _1 . Above 
shock speeds of about 35 km s _1 the highest sputtering rate per 
unit time per unit volume of elemental silicon is due to collisions 
with H2O. Table [3] displays the fraction of the initial core elemen- 
tal silicon that is sputtered. It also gives the enhancement in the 
downstream gas phase SiO fractional abundance due to sputtering. 

Similarly to the results for mantle sputtering, the core sputter- 
ing of the small grains is independent of nno (as their Hall parame- 
ter is much larger than unity). The large grains are more efficiently 
sputtered in the low-density regime but contribute only a negligi- 
ble amount to the SiO gas-phase abundance. For a range of shock 
speeds, the SiO producti on rises with i ncreasi ng speed with val- 
ues similar as in Fig. 2 of lGusdorf etafl j2008j) . However, the SiO 
abundance reaches an upper limit as C-type shocks do not exist 
for shock speeds above a certain value. For higher shock speeds, 
shocks are of J-type. For tiho = 10 4 cm -3 the transition occurs 
for a shock speed between 50 and 55 km s , while a J-type shock 
occurs for a minimum shock speed between 35 and 40 km s _1 for 
tiho = 10 6 cm~ 3 . At these speeds the grain cores are not com- 
pletely destroyed, i.e. only a few percent of the core is removed in 
the high-density case and about a third of the core is destroyed in 
the low-density regime. Hence, the maximum SiO fractional abun- 
dance differs by an order of magnitude between the two density 
regimes. 



3.4 Angular dependence 

Up to now we have discussed the g rain sputtering for an oblique 
shock for which 9 = 45°. However. I Van Loo et al.l J2009I ) showed 
that the shock properties, such as the shock width, change as func- 
tion of 9. Therefore, variations, with angle, of the H2O and SiO 
abundances due to sputtering is expected. 

Figure [2] presents results for 16 km s _1 shocks propagating 
into media with n H o = 10 4 cm~ 3 at angles of 30°, 45° and 60° 
with respect to the upstream magnetic field. The fraction of the 
ice that is sputtered from each type of grain, rather than the gas 
phase abundance of water, is shown because much of the gas phase 
water exists due to the conversion of upstream gas phase atomic 
oxygen to water. The results in Fig. [2] indicate that the H2O abun- 
dance is lowest in the shock with the largest width. It is obvious 
that the shock width is irrelevant for the level of sputtering as the 
sputtering depends mainly on the neutral-grain drift speeds. As the 
bottom row of Fig.[2]shows, the peak small and large grain-neutral 
drift speeds are smallest at large angles; they are 10.5 km s _1 at 
30°, 10 km s _1 at 45° and 9.5 km s _1 at 60°. This behaviour is 
consistent with the results shown in Fig. 2 o f lVan" Loo et al. U2009T) 
who found that the y and z components of the velocities of charged 
species are greater for smaller angles. Furthermore, the length scale 
over which the drift speeds are higher than 6 km s~ 4 , i.e. the mini- 
mum drift speed for which mantle sputtering is efficient, are similar 
for all angles. The angular variations in the fraction of water that 
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Table 3. Fraction of core elemental silicon that is sputtered and gas phase SiO fractional abundance enhancement due to sputtering. 
Results are given for large and small grains and for both combined. 
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2.20(-6) 


45 
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1.15(-6) 


3.04(-l) 
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9.28(-2) 


5.40(-6) 


20 


1.00(6) 


0.00 




0.00 




0.00 




25 


1.00(6) 


0.00 




7.37(-6) 


1.03(-10) 


1.77(-6) 


1.03(-10) 


30 


1.00(6) 


0.00 




3.22(-3) 


4.50(-8) 


7.73(-4) 


4.50(-8) 


35 


1.00(6) 


9.76(-4) 


4.31 (-8) 


3.26(-2) 


4.56(-7) 


8.57(-3) 


4.99(-7) 




Distance (10 15 cm) Distance (10 15 cm) Distance (10 15 cm) 



Figure 4. Gas phase fractional abundance of SiO for a shock speed of 25 km s — 1 , njjo = 10 6 cm -3 and 9 of 30° (left column), 45° (middle column) and 60° 
(right column). The top row gives the ion (black, dotted), electron (red, dashed) and neutral (blue, solid) temperatures, as well as the neutral density normalised 
to its upstream value (green, dash-dotted). The middle row shows the gas phase fractional abundances of SiO. The bottom row displays the speed of small 
(red, dotted) and large (blue, dashed) grains relative to neutrals in km s — 1 . 
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is sputtered are significant and range from 0.4 at 60° to 0.6 at 30° 
Such variations disappear with increasing shock speed because for 
high enough speeds vaporisation of the grain mantle is complete 
for all angles. 

Figure [4] presents results for 25 km s~ shocks propagating 
into media with n H o = 10 6 cm -3 at angles of 30°, 45° and 60° with 
respect to the upstream medium. The results in Fig.[4]show that the 
SiO abundance too is lowest in the shock with the largest width. 
Again these differences can be attributed to the peak small grain- 
neutral drift speeds which are smallest at large angles; they are 18.2 
km s" 1 at 30°, 18 km s" 1 at 45° and 17 km s" 1 at 60°. Contrary 
to the H2O results, these modest differences in drift speeds have 
large effects on the amount of elemental silicon that is sputtered, 
i.e. there is a difference of two orders of magnitude between the 
SiO abundance at 30° and 60° . However, note that the shock speed 
is close to the threshold speed for core sputtering. Then small vari- 
ations in the drift speed can mean the difference between core sput- 
tering or no sputtering at all. For higher shock speeds the angular 
differences in SiO abundance are less pronounced, but the trend of 
lower abundances for larger 8 remains. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we present the first self-consistent multifluid MHD 
simulations of grain mantle and core sputtering in oblique C-type 
shocks propagating through dusty media. The small grains are more 
tightly coupled to the magnetic field than the large grains with a 
strong coupling in both low- and high-density regimes. This max- 
imises the drift speeds between the small grains and the neutral 
sputtering projectiles resulting in more efficient sputtering of the 
small grains. Mantle sputtering starts at shock speeds of about 
10 km s" 1 , and the mantles are completely vapourised for shock 
speeds around 20-25 km s" 1 . The sputtering of the icy mantles re- 
turns water back to the gas phase. However, neutral-neutral chem- 
istry dominates the production of gas phase water by an order of 
magnitude. 

Once the mantles are completely eroded, core sputtering com- 
mences and releases silicon to form SiO. The abundance of SiO 
in the gas phase never saturates as even at the critical speed for 
the transition from C- to J-type shocks not all of a grain core is 
destroyed. While the critical speed is 50-55 km s" 1 for tiho = 
10 4 cm -3 , it is only 40 km s" 1 for who = 10 6 cm" 3 . For a given 
shock speed the efficiency of small grain sputtering is not sensitive 
to the density. Consequently, the density dependence of the C- to J- 
type transition means that higher SiO abundances can be produced 
by shocks in the low-density regime. 

An important result is that the sputtering yields depend on 
the angle between the magnetic field and the shock normal. Es- 
pecially, at the threshold drift speeds for mantle or core sputter- 
ing, small angular variations can induce large changes in abun- 
dances. For example, for shocks propagating at 25 km s~ 1 in media 
with number densities of 10 6 cm -3 , the SiO abundance differs by 
two orders of magnitude for 8 — 30° and 60°. This has impor- 
tant consequences for the interpretation of observations. For exam- 
ple, discrete emission features in SiO observations have been as- 
crib ed to different shocks bein g present within the telescope beam 
(e.g. ljimenez-Serra et al.|[2009h . However, such features can be ex- 
plained by a bow shock interacting with an inhomogeneity in the 
medium. The flanks of the bow shock have lower speeds than the 
head and also make different angles with the magnetic field. In an 
inhomogeneous medium some regions of the bow-shock are then 



highlighted resulting in discrete SiO emission features. Such an as- 
sumption, however, needs to be further examined by extending the 
planar shock model to multiple dimensions. 

Direct comparisons of the results with those of other re- 
searc hers is complicated as the only previous work dCaselli et al.l 
1 19971) on sputtering in oblique shocks is not based on self- 
consistent shock models. Also, the treatme nt of core sputteri ng in 
that work differed from that developed by I May et all fcOOOl) and 
used in subsequent work. However, the results given in Table[3]for 
the fraction of elemental silic on sputtered from s mall grains are 
in harmony with the results o f lCaselli et all jl997t) for a shock for 
which 9 — 45°. Other self-consistent studies concern only per- 
pendicular shocks (8 — 90°) and, thus, should find similar, al- 
though lower, values for shock speed exceeding the threshold value 
for sputtering. iGusdorf et al.l d2008al) . who assumed a somewhat 
larger upstream magnetic field strength and different grain popula- 
tion, found that the fraction of silicon that is sputtered is about 0.04 
fornHo = 10 4 cm s _1 and a shock speed 40 km s ~ 1 . iGuillet et al.l 
J201 lj) who adopt a kinetic method for the dust grain dynamics in 
fa vour of a fluid approac h, find a value of 0.02 for the same model 
as lGusdorf et al.l J2008ah . The corresponding result given in Table 
Sis 0.04. 

The SiO abundances found in our simulations only repre- 
sent upper limits. Observations indicate that silicate-based grains 
only account for about 75% of the interstellar dust mass with the 
rest o f the mass stored in carbonaceous grains or PAHs jDrainel 
l201ll) . Furthermore, the conversion of atomic silicon into SiO is 
not fully efficient above densities of 10 s cm" 3 which possibly re- 
duces the_Si_0_abundance due to sputtering by an order of magni- 
tude l lGusdorf et ai]|2008bl) 

Unlike IGusdorf et ai] j2008j) and IGuillet etail feOllh . who 
included grains 40 times smaller than the smallest that we have 
assumed to be present, we have not considered the destruction of 
grain material due to grain-grain collisions. Our choice to focus 
solely on sputtering first was based in part on a desire to under- 
stand the effects of obliqueness on sputtering before looking at the 
complicated way in which obliqueness, sputtering and grain-grain 
collisions all act together to dest roy grains. Our deci sion to do so 
was informed by the estimates of ICaselli et alj J 19971) that showed 
for grain populations of the sizes that we selected, grain-grain col- 
lisions would contribute significantly to mantle removal only at 
shock speeds near the threshold even for high pre-shock densities. 
Our choice has allowed us to see that obliqueness affects sputtering 
near the threshold s hock speed significa ntly. However, cons i dering 
the findings of both lCaselli et al.l ill 9971) and o flGuillet et ail feoill) 
for the return of silicon to the gas phase, as well as the results given 
in Figure [4] we conclude that the future inclusion of both grain- 
grain collisions and obliqueness in shock models is necessary to 
infer the threshold shock speed for the injection of significant quan- 
tities of silicon into the gas phase. 
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APPENDIX A: PARTICLE CURRENT ONTO A GRAIN 



lHavnes etaiUl987h derived general expressions for the current of 
charged particles (with charge Zj) onto a dust grain carrying a 
charge of Z g . 



For Z g Zj «: 0, 



z *£-{( 1+2 [ v i) 



x erf ( m. \ + ^HmM 



exp 



and, for Z q Z, ^ 0, 



Tj g — Tia rijZje 



j fi 

2v la I 2 



1 + 2 



ZZjZgS, 



f P + yW + erf — -v/W 



+ 



+ v/^J^ 9 C exp 



- JZjZaJ 



where a is the grain radius, rij the number density of the charged 
fluid, Vj g the relative drift velocity between the grain and the 
charged fluid, Cj the sound speed of the charged fluid and £ = 
e 2 jaksT. The charged fluid can be the ion or electron fluid. 
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